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Abstract 

We develop statistical mechanics and thermodynamics of Bose and Fermi systems 
in relativistic harmonic oscillator (RHO) confining potential, which may be applicable 
in quark gluon plasma (QGP), astrophysics, Bose-Einstein condensation (BEC) etc. 
Detailed study of QGP system is carried out and compared with lattice results. Further, 
as an application, our equation of state (EoS) of QGP is used to study compact stars 
like quark star. 
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1 Introduction 



There are many systems like QGP, quark stars, BEC etc. where particles are confined by 
harmonic oscillator type potentials and are at thermodynamic equilibrium. So we need to 
study statistical mechanics of such a system. Taking QGP, a new form of matter made of 
quarks and gluons, as an example we discuss this problem. It is generally accepted that 
hadrons are confined state of quarks and gluons and at high temperature and density, we get 
what is called QGP [1]. In relativistic heavy ion collisions (RHICs), quark stars etc., QGP 
consists of quarks and gluons confined in a volume large compared to the size of hadrons. In 
the case of early universe, just before hadronization, again we have QGP. In all above cases, 
EoS of QGP is needed to explain their evolution and study various properties. Since quantum 
chromodynamics (QCD) is nonlinear, it is difficult to solve analytically at low energy and 
only at high energy one can do perturbative calculations. That is why one goes to numerical 
simulations of QCD on lattice and EoS of QGP without and with quark flavors has been 
derived [2]. Many phenomenological models like hard thermal loop (HTL) [3], quasi-particle 
QGP (qQGP) [4, 5], strongly interacting QGP (sQGP) [6], strongly coupled QGP (SCQGP) 
[7] and various confinement models [8, 5, 9] were proposed to explain lattice results and 
there is no satisfactory model. Models which fits well with lattice results, generally, involve 
large number of parameters [4, 5] and models with detailed theory like HTL [3] don't fit 
lattice results well, so on. However, recently we found that SCQGP model [7] as well as a 
new qQGP model [10] fits surprisingly well with lattice results with minimum number of 
parameters compared to other models. Both models are purely phenomenological and based 
on the collective properties of plasma. 

Here we approach the problem from the known fact that RHO model of hadrons quite 
successfully explains many properties of hadrons [11]. QGP formed in RHICs and compact 
stars may be viewed as a large hadron with large number of quarks and gluons, confined in 
a RHO confining potential. Hence, a study of statistical mechanics and thermodynamics of 
such a system will give EoS of QGP. It is just a generalization of widely studied Bag Model of 
QGP. Instead of sharp bag wall, here we have a smooth confining potential. As a first step, 
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here we neglect effects like Color-Coulomb interactions, hadron formations etc. and see how 
the EoS of QGP differs from ideal EoS. Earlier, similar calculations has been done [9], but 
without proper counting of degeneracy. They have taken only accidental degeneracy where 
as here we have taken full degeneracy arising from coalese of continum levels to oscillator 
levels as a result of RHO interaction, just like in Landau's theory of diamagnetism [12]. 

2 EoS of QGP in RHO Model : 

From the earlier study of hadrons in RHO model, we have the single particle energy levels 
[11], 

el = C 2 g (2n+l), (1) 
for gluons and 



e 2 n = m 2 q + C^e n + m q (2n + l), (2) 

for quarks. C g and C q are spring constants of gluon and quark RHO confining potentials 
respectively, n = 1,2,... is the oscillator quantum numbers and m q is the mass of quark. 
Let us first consider gluon plasma without quarks which is a Bose system. Following the 
standard procedure of statistical mechanics, pressure or logarithm of partition function is 
given by, 

^ = -E^/m(l-e-^), (3) 

n 

where g n is the full degeneracy arising from the "coalescing together" of an almost continuous 
set of the zero-potential levels and is given by, 

vr 3 

S n = -^((2n + 3) 3 / 2 -(2n+l) 3 / 2 ). (4) 

V is the volume and (3 is the inverse temperature, gi = 16 is the internal degrees of freedom 
of gluons. Since gluons are bosons, the fugacity z — 1. For fermions like quarks z ^ 1 and 
pressure is given by, 

^ = £ff»0/M 1 + * e " /fc ")> (5) 



where gj is the internal degrees of quarks. Above expressions for pressure are exact and one 
need to do infinite sum which may be done numerically. However, for high temperature, we 
may use Euler-Maclaurin series to get an appproximate expression as, after some algebra, 



■ CgK i 3 cy ^ i , C g L 

TAT ^ T ' 



+^(3V3-l)Ml-e-*) (6) 

for gluon plasma, in terms of modified Bessel functions K v (x). But we know that at ex- 
tremely high temperature, confinement effects may be negligible and we expect to get EoS 
of relativistic free gas which we indeed get as shown below. As T — > oo, 

90 An 2 T 2 { ' 

^ 9 -^T^P s . 
90 

Note that in our quantum calculation with proper choice of g n , Eq. (4), as expected, we get 
T 4 , rather than T 7 as in [9]. Further, we also know that as T — > oo, we approach classical 
limit and from our earlier classical calculations of this problem, we indeed got T 4 law [13]. 
Next let us compare our EoS with the results of lattice simulation of QGP or lattice gauge 
theory (LGT). 

3 Comparison with LGT: 

We can see, from Eq. (7), that our EoS of QGP is temperature dependent modification of 
normal T 4 law. Similar modifications is seen in LGT calculations and only at T — > oo, one 
expects T 4 law, which is termed as non-ideal effect. As an example, let us consider pure 
gauge QGP or gluon plasma and the pressure is, Eq. (3), 
P 

T ~~ 6tt 2 
which may written as 

p - y s = ¥ E((2n + 3)3/2 " (2n + 1)3/2) ln(1 " e ~ f ^ ' (9) 



2((2n + 3) 3 / 2 - (2n + if' 2 ) ln(l - e-^v^+T) ? (8) 



where P s is the Stefan-Boltzman limit pressure and t = T/T c . T c is the critical temperature 
of the transition of hadrons to QGP and the parameter a is defined as, a = C g /T c . Above 
equation may be compared with LGT results by varying the parameter a as discussed in 
section V. Knowing P or partition function, all other thermodynamic quantities like energy 
density, entropy, sound speed etc. may be derived. Further extension of our calculations 
for QGP with finite number of quark flavors, without or with finite chemical potential, is 
straight forward [14]. 

4 Quark star: 

As another example, let us consider what are called quark stars [15]. We know that, at the 
last stages of a star, it undergoes various changes like white-dwarf star, neutron star and 
then, probably, quark star. Sometimes more complex compact star with different layers of 
matter with quark matter inside and then mixture of quarks and hadrons, neutrons, so on. 
Here we take a very simple compact star, just cold quark star, where the degenerate pressure 
of quarks in RHO model balances the gravitational collapse. We estimate the limiting mass 
for stability and mass-radius relation, using Newtonian description. A detailed study using 
Tolman-Oppenheimer-Volkoff (TOV) equation with complex compact star in our model is a 
separate topic. With these simplifications, we get the degenerate pressure from the energy 
of the system, 



u " n=l 

where we have approximated the single particle energy level of quarks, Eq. (2), for simplicity, 

as [9] 



Here gf is the internal degrees of quarks and m q = rip is the quantum number corre- 
sponding to Fermi level, given by the relation, 




(10) 



e 2 n = m 2 q + C 2 q {2n+l). 



(11) 



N 




£((2n + 3) 3 / 2 -(2n + l) 3 / 2 ), 



(12) 



6tt 2 



n=l 
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where N is the total number of quarks and in the case of cold star approximation, quarks 
fill the single particle RHO levels, according to Pauli's exclusion principle, upto the level %. 
Following the same procedure as that of white-dwarf star calculations we may relate the n F 
to the radius, R, of the star as 

R = f 9 * M - J , (13) 

where M is the mass of star. Degenerate pressure is obtained from E , 
9E 



dV 

Cj9f 
6tt 2 



1 



(2n F + 3) 3 / 2 - 3^ d 



3{2n F + 3)V2 dn F _ 



np 

^((2n+3) 3/2 -(2n+l) 3/2 ) + 2n + 1 .(14) 



71=1 



To proceed further, we approximate the summation over n by integration using Euler- 
Maclaurin series. The error introduced by this approximation is small, which may be checked 
numerically by evaluating summation and Euler-Maclaurin formula. This degenerate pres- 
sure balances the gravitational pressure inward, which we take, the Newtonian formula, 

ai M 2 

t^w < 15 > 

where 7 is the gravitational constant and a is a geometrical factor. After some lengthy 
calculations for large n F , we get 



(i? 2 + l)(2-i? 2 ) 2 = M 4/3 , 
where R and M are related to radius, R, and mass, M, of the star respectively as, 



R — R C„ 



■ 2 9f ra g i /3 
- 9ttM ' 



and 



M 



M 



[a j ml) 



2yi/2 / 9f ^ 
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(16) 



(17) 



(18) 



It is interesting to compare above mass-radius relation with familiar white-dwarf star mass- 
radius relation, proceeding with similar normalization of mass and radius, 



( J R 2 + l)(2/3-i? 2 ) 2 = M 4 / 3 . 



(19) 



We see that the only change in our case is that the factor 2/3 is replaced by 2 in mass- 
radius relation and C q plays the role of m e , the electron mass. From the Eq. (16), we can 



where m p is the mass of proton. It is generally assumed [15] that quark star contains up, 
down and strange-antistrange quarks and hence ~ (2x3x2 + 1x2x3x2) = 24. The 
quark mass m q in RHO is a constituent mass which ranges from 300 MeV to 600 MeV for 
up quark to strange quark and on an average we may approximate m q ss m p /2. Therefore, 
as a rough estimate, from Eq. (20), limiting mass of quark star may be of the order of 4 
times the solar mass, where as the Mq is of the order of solar mass. Further numerical 
calculation of Eq. (16), without the approximation of rip — > oo which is used to get the Eq. 
(16), we obtained a plot of mass-radius relation as shown Fig. 2. 



In Fig. 1, we plotted our results for pressure in terms of p = P/P s as function of t = T/T c 
along with LGT results. By varying the parameter a, so as to get a similar behaviour as 
LGT results as function of t, we get a = 2.3. However, the curve lies slightly above the 
lattice data. We may correct our pressure by subtracting a constant A such that both the 
curve matches at large t values and we get good fit except near to T = T c . The value 
of A here is 0.1. It is not unreasonable to subtract by A because we have neglected the 
mutual interactions between particles, mainly, attractive Color-Coulombic interactions. We 
know that attractive mutual interactions among particles reduces the pressure. From the 
expression for a, we may estimate T c = C g /a 200 MeV" for C g = 2.35 jmT 1 [9]. Once 
pressure or partition function is obtained all other thermodynamic functions like e, c 2 s etc. 
may be easily evaluated. Similar comparison of our results and LGT results may be carried 
for other systems with flavored QGP [14]. 



estimate the limiting mass, M® s , similar to the Chandrasekhar's limit (M^ 7 ) in the case of 
white-dwarf star, above which the star will collapse and we get, 




(20) 



5 Results : 
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In Fig. 2, we plotted mass-radius relations of quark star. Mass is expressed in terms of 
solar mass and radius in kilometers (km) by solving, numerically, Eq. (14) after equating 
to the gravitational pressure. We have taken m q = 2.5 which is equivalent to taking C q ~ 
lfm- 1 [9]. 

6 Conclusions : 

We developed statistical mechanics and thermodynamics of a system of quantum particles in 
RHO confining potential and applied to QGP. We compared the results with LGT and con- 
clude that the overall QCD confining effect, modelled by RHO, explains the LGT data with 
small correction to take into account of the attractive mutual interactions among particles. 
We applied our model to study quark star and obtained the limiting mass and mass-radius 
relation and found to agree of the order of magnitude. To get an accurate result one need to 
study complex compact stars using TOV equations. However, Boson or Fermion system in 
RHO confining potential is a new general problem which may be applicable to many fields 
like BEC, RHICs, astrophysics and so on. 
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Figure 1: Plots of P/P s as a function of T/T c from our model and lattice results (symbols) 
for pure guage. 

Figure 2: Mass-radius plot of quark star, mass in terms of solar mass and radius in kilometers. 
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